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Abstract The magnetospheres of neutron stars provide a valuable testing ground 
for as-yet unverified theoretical predictions of quantum electrodynamics (QED) in 
strong electromagnetic fields. Exhibiting magnetic field strengths well in excess of 
a TeraGauss, such compact astrophysical environments permit the action of ex- 
otic mechanisms that are forbidden by symmetries in field-free regions. Foremost 
Ph among these processes are single-photon pair creation, where a photon converts to an 

electron-positron pair, and magnetic photon splitting, where a single photon divides 
W~) into two of lesser energy via the coupling to the external field. The pair conversion 

, , process is exponentially small in weak fields, and provides the leading order con- 

tribution to vacuum polarization. In contrast, photon splitting possesses no energy 
threshold and can operate in kinematic regimes where the lower order pair conversion 
is energetically forbidden. This paper outlines some of the key physical aspects of 
these processes, and highlights their manifestation in neutron star magnetospheres. 
i Anticipated observational signatures include profound absorption turnovers in pulsar 

spectra at gamma-ray wavelengths. The shapes of these turnovers provide diagnos- 
tics on the possible action of pair creation and the geometrical locale of the photon 
emission region. There is real potential for the first confirmation of strong field QED 
with the new GLAST mission, to be launched by NASA in 2008. Suppression of 
pair creation by photon splitting and its implications for pulsars is also discussed. 
Keywords : Quantum Electrodynamics, Magnetic Fields, Neutron Stars, Pulsars. 

OO 

1 Introduction 
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^ Quantum Electrodynamics (QED) is one of the most robust theories of physics. It 

has been extensively and almost exhaustively tested in terrestrial laboratory experi- 
ments, as have its constituent, ingredient disciplines, special relativity and quantum 
mechanics. The thoroughness of its investigation has forged an extreme confidence 
in the physics community in the viability of its predictions, and its broad applica- 
bility in Nature. This has naturally extended to astrophysical settings, where both 
the more conventional and even the more exotic processes from QED are applied 
in models of different sources of radiation. Foremost among these are Compton 
scattering — > e7 and two-photon pair creation, 77 — > e + e~ , which have been 
precisely calculated and tested since the 1930s as processes that are probable in the 



cosmos due to their comparatively high cross sections. Because they most gener- 
ally apply to relativistic environs where X-rays and gamma-rays abound, they are 
ubiquitous in models of compact, non- magnetic astrophysical systems, particularly 
those involving black holes, both galactic (typically with masses M of a few solar 
masses M & ), and extragalactic and supermassive (with M ~ 10 6 — 10 9 M Q ). 

Neutron stars, compact remnants of aged stars more massive than the sun that 
have undergone a supernova event, provide a very different environment for physical 
processes. Supported only by neutron degeneracy pressure against the incredibly 
compressing pull of gravity, they stably exist with radii i?NS of the order of 10 km. 
During the supernova core implosion, it is believed that charge currents in their 
outer layers persist and strengthen so that their magnetic fields intensify, perhaps 
exceeding magnetic flux conservation. They can then possess fields in the range of 
10 10 - 10 15 Gauss (i.e., 10 6 - 10 11 Tesla; see, e.g. [H |2]) that are far beyond the 
realm of terrestrial experiments. The prime observational manifestations of neutron 
stars are pulsars, rapidly pulsating radio sources that were discovered in 1967 [3]. 
Now seen in other wavebands such as the optical, X-rays and gamma-rays, their 
interpretation as rapidly rotating neutron stars [I] that slowly spin down to longer 
pulse periods has provided powerful evidence for the existence of strong B fields. 
This evidence is derived from a rotating, tilted dipole model for the field structure 
[6], the rotational/electromagnetic torque on which leads to inferences of surface 
fields in the range B ~ 10 9 — 10 15 Gauss. Such a dipole picture is the simplest 
choice for the field geometry, being the leading order moment for any magnetic 
field configuration. The rotation and magnetic axes are generally not aligned, just 
as in the cases of the sun and the Earth, and so the time-varying B fields induce 
electric fields E of substantial magnitude above the stellar surface. These E fields 
are quenched [6] on relatively small lengthscales by charges that are prolifically 
created via the magnetic pair production process that is detailed in this paper. 
Furthermore, note that additional support for the existence of such intense B fields 
comes from spectroscopic observations of putative cyclotron emission or absorption 
lines in the 20-150 keV band of accreting X-ray binary pulsars [7]. A discussion of 
these pertinent astrophysical elements can be found in [2]. 

These strong field environs provide an opportunity for more exotic and less fa- 
miliar predictions of QED to come to the fore. The B field profoundly changes the 
physical nature of quantum interactions, yielding a strong bias towards momentum 
exchange along the field. It quantizes the electron states in a direction perpendicular 
to the field so that momentum in that direction is no longer conserved; a product 
of the lack of translational invariance of the states orthogonal to B. Accordingly, 
B = symmetries are broken, and interactions that are forbidden in field-free en- 
virons become possible, and even probable in neutron star magnetospheres. The 
most likely among these are magnetic pair creation 7 — > e + e~ and photon splitting 
7 — > 77 . These are more recent predictions of QED, dating from the 1950s and 
1960s, partly due to their inherent mathematical complexity (spawned by quanti- 



zation in a cylindrical geometry), and partly due to the contemporaneous rising 
interest in their applicability to cosmic systems, e.g. neutron stars. The fact that 
these, as yet, untested predictions of QED are widely invoked in pulsar contexts 
is a true testament to physicists' confidence in the theory of relativistic quantum 
mechanics. These two processes form the focus of this paper, though it is noted that 
strong-field quantum effects also have profound influences on Compton scattering 
[8j [9] , classical synchrotron radiation [101 E] , and other radiation processes. 

2 Magnetic Pair Creation 

One-photon pair production 7 — > e + e~ in strong magnetic fields is a first-order 
QED process with one vertex that is quite familiar to pulsar theorists, having been 
invoked in polar cap models to explain the photon emission of both radio pulsars 
(e.g. [5j [6j) and the handful of known gamma-ray pulsars (e.g. [12]). It is for- 
bidden in field-free regions due to the imposition of four-momentum conservation, 
but takes place in an external magnetic field, which can absorb momentum perpen- 
dicular to B (momentum along the field is still conserved, as is the energy). The 
first analytic computations of its rate R pp [131 EI] (hereafter the superscript pp de- 
notes pair production) indicated a rapid rise with increasing photon energy e 7 and 
magnetic field strength, becoming significant for 7 -rays above the pair threshold, 
uj = 6 y /(m e c 2 ) = 2/ sin#kB > an d for fields approaching the quantum critical field 



Here #kB is the angle the photon momentum vector k makes with the magnetic field 
B. B cr represents the field at which the cyclotron energy equals m e c 2 , and defines 
the field scale at which the impact of the external B field on quantum processes 
becomes significant. It is appropriate to identify now the dimensionless scaling 
conventions to be adopted throughout this paper that are of common usage. All 
photon energies, u , will be scaled in terms of the electron rest mass energy m e c 2 , 
and all field strengths will hereafter be expressed in terms of B cr , i.e. B = 1 denotes 
a field of 4.413 x 10 13 Gauss. We note in passing that pair creation 7 — > e + e~ 
in the induced electric fields is also possible in rotating pulsar magnetospheres, 
however this is reduced by the order of v TOt /c <C 1 near the neutron star surface 
due to corotation speeds v TOt = 2irR m /P usually being relatively small ( i? NS is the 
neutron star radius and P is the pulsar period). 

The magnetic pair creation rate is resonant at the thresholds for each combina- 
tion of produced pair states, due to the available momentum parallel to the field 
approaching zero in the frame where = tt/2 . A large number of integrable (over 
photon energies) resonances result, producing a characteristic sawtooth structure 
[T5| ITS] that is displayed in Figure [l] below. In astrophysical contexts, a range of 
field strengths and initial photon energies u are always sampled, and the resulting 
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convolution smears out the sawtooth appearance into a continuum. Fully general 
expositions of the 7 — ► e + e~ rate in uniform B fields can be found in [T(H [15]. Here 
we highlight the spectral structure near absolute pair threshold uj = 2/sin6* kB , 
a case most germane to pulsar applications since generally photons are produced 
in beams almost along the local field lines in the dipole geometry, so that pair 
creation threshold is crossed from below during magnetospheric propagation of pho- 
tons. This implies that calculations in the domain uj ^ 2/ sin^B are typically of 
greater practical relevance, especially for near-critical or supercritical fields B ^ 0.1 
[TT] . The pertinent rates, exhibited in [T7j, can be written as follows. Let (j,k) 
denotes the Landau level quantum numbers of the produced pairs, which have ener- 
gies Jl +pj}~ + 2jB and + p 2 k + 2kB . The parallel momenta pjk of the pairs 
in the frame of reference for photons moving perpendicular to the field are given by 



Pjk 



^sin 2 fl kB -l-(j + fc)i?+ ^-^ X2 12 



uj sin #kB 



(2) 



with the solution of pjk = defining the host of resonant energies uj . Throughout, 
electron energies are scaled by m e c 2 and their momenta are in terms of m e c . 

Near threshold, only a few pair states are kinematically accessible. The rates are 
dependent on the polarization of the incoming photon. Following common practice, 
here we adopt the convention that photon linear polarizations are such that || refers 
to the state with the photon's electric field vector parallel to the plane containing 
the magnetic field and the photon's momentum vector, while _L denotes when the 
photon's electric field vector is normal to this plane. Let A = h/m e c be the Compton 
wavelength of the electron over 2tt , and a { = e 2 /he be the fine structure constant. 
The exact, polarization-dependent, pair production rate [IS], including only the 
(j = 0, k = 0) pair state for || polarization is: 

for £ = u 2 sin 2 O^b/^B] , while the sum of the (j = 0, k — 1) and (j — 1, k — 0) 
states contributes just above threshold for the _L polarization: 



- r W ( E ° E > + 1 + A) exp(-0 , a, > - (4) 

The pertinent dimensionless energies of the produced pairs are 



E = ^l+p 2 01 , Ei = y 1 + P01 + 2-B • (5) 

Observe that the absolute threshold for _L photons is higher than that for || photons, 
an effect that finds its origin in the spin-dependence of energies in the final pair 
states. Indeed, the rate for the _L state is always lower than that for the || state. 
This property for the leading order photon absorption process, when connected via 



the optical theorem to the consideration of dispersion in the magnetized vacuum, 
clearly indicates that such a vacuum is birefringent. 

For photon energies well above threshold, the number of accessible pair states 
becomes very large, much larger than just the few states considered above (see 
Figure[l]for an illustration of the sawtooth resonant structure of the rates). Then the 
resonances at pjk = tend to blend into one another, an effect that is enhanced in 
pulsar magnetospheres by sampling significant ranges of photon energies, angles and 
field strengths that smear out the cusp-like structure in forming an average reaction 
rate. In this regime, the dependence on u can be averaged over ranges larger than 
the typical separation of the resonances, and one can use more convenient asymptotic 
expressions for the polarization-dependent pair creation rates [HI HH EE] : 

i?J p « 2Rf « l^l^sinflkB expj-J^}, X = u Bern fa < 1 , (6) 

where x is the critical asymptotic expansion parameter. These results essentially 
derive from Schwinger-type mathematical formulations, as do alternative results 
(e.g. [18j [19]) for the x ^ 1 regime that are generally of less practical interest. 
The differences between the asymptotic result in Eq. ^ and the exact rates be- 
come profound near threshold [T3], and the asymptotic analyses can be improved 
considerably by treating mildly-relativistic regimes for the produced pairs [TTl [To] . 

As mentioned above, in polar cap pulsar models [5], El [T2j , high energy radiation is 
usually emitted at very small angles to the magnetic field, well below pair threshold 
[T7] . The 7 -ray photons will convert into pairs only after traveling a distance s that 
is a fraction of the field line radius of curvature p c (which exceeds the neutron star 
radius _R NS ), so that sin# kB ~ s/p c . From Eq. ([6]), the pair production rate will be 
vanishingly small until the argument of the exponential approaches unity, i.e., when 
cu-Bsin^kB ~ 0.2. Consequently, pair production will occur well above threshold 
when B ^0.1 and the asymptotic expression will be valid; from Figure [T] one can 
deduce that the attenuation length will then be much less than i? NS . In contrast, 
when B <; 0.1 , pair creation will occur at or near threshold, and Eqs. ^ and Q 
become more appropriate. In this domain, an additional sublety arises, since another 
mode of pair creation exists, namely the formation of pairs in a bound state, i.e. 
positronium. This has been proposed as an effective competitor to the production 
of free pairs [201 IZTJ [22] because the binding energy lowers the threshold slightly 
(< 1%) below the value for production of free pairs. The relevance of positronium 
formation to pulsars is potentially great if the bound state is stable for considerable 
times. Positronium is subject to destruction by three main mechanisms: free decay, 
electric field ionization, and photo- ionization. As discussed in [T7] , due to conflicting 
calculations, it is presently unclear whether or not positronium is stable to free decay 
and photo-ionization; assessment of this issue motivates future work. 
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Fig. 1: Attenuation coefficients or inverse attenuation lengths (i.e., rates i?^ p and 
-Rj^im divided by c), for pair creation (above the threshold u = 2) and photon 
splitting for the polarization mode -L— (only for uj < 2), as functions of the 
incident photon energy uj , for field strengths B , as labelled. Photons are assumed 
to propagate orthogonally to the field lines ( #kB = 90° ). The pair creation rates are 
computed using [151 123], and only for _L photons; at to > 4 the dashed B — 0.1 
curve is the asymptotic form in Eq. (|6]). Also, the dashed splitting curve labelled 
B = 10 3 represents the asymptotic high-field "saturation" result. The horizontal 
heavy dashed line labelled 1/Rns marks the approximate minimum effective attenu- 
ation coefficient for which these processes become prolific in pulsar magnetospheres. 

3 Magnetic Photon Splitting 

The most important competitor to pair creation 7 — > e + e~ at high magnetic field 
strengths as a mechanism for attenuating photons in pulsar magnetospheres is mag- 
netic photon splitting 7 — > 77 . The relevance of this process to neutron star envi- 
ronments has been emphasized in various papers [Ml [25], EH] • Splitting is a third- 
order QED process with a triangular Feynman diagram, and therefore is weaker 
than magnetic pair creation by the order of aj , yet without the constraint of a 
threshold energy, an important distinction. Though splitting is permitted by energy 
and momentum conservation, when B = it is forbidden by a charge conjugation 
symmetry of QED known as Furry's theorem (e.g. [2S]), which states that ring dia- 
grams that have an odd number of vertices with only external photon lines generate 



interaction matrix elements that are identically zero. The presence of an external 
field breaks this symmetry. The splitting of photons is therefore a purely quantum 
effect, and has appreciable reaction rates only when B ^ 1 . 

Magnetic splitting 7 — > 77 is a relatively recent prediction of QED. A decade 
of controversy followed the earliest calculations in the 1960s before the first correct 
evaluations of its rate [2U E7J [28] were performed, via an effective Langrangian tech- 
nique. These works focused on asymptotic forms in the limit of photon energies well 
below pair creation threshold, which varied as B 6 when Bel. Their rate deter- 
minations neglected dispersion in the birefringent, magnetized vacuum, so that all 
photon momenta were assumed collinear. The early controversy was fueled by the 
inherent difficulties in calculating the rates of this third order process by standard 
QED techniques, an issue that re-emerged in the 1990s (see [23] for a discussion). 
For weak-dispersion regimes, fully general rates in Schwinger-type formalisms were 
derived by various groups [2U ESI EDI EI]- Even more general S-matrix, Landau 
representation calculations that include information on the pair resonances are de- 
rived by [32] and [331 El] (note that issues concerning the erroneous numerics in 
[32] are addressed in [23]). These resonances are a significant factor in determining 
the splitting rate near and above pair threshold, since then the intermediate pair 
states go "on-shell" and the splitting process effectively becomes first-order in a { , 
tantamount to a pair conversion followed by cyclotron decays. 

It is practical to restrict considerations of splitting to regimes of weak disper- 
sion, where manageable expressions for its rates are obtainable. These are still 
complicated triple integrations in Schwinger formalisms [2H [291 EOl El] or triple 
summations in S-matrix calculations |34| . Further specialization to either low mag- 
netic fields ( B <C B cr ) or low photon energies (w<2) therefore proves expedient, 
and palatable results for splitting rates were first obtained in such regimes. A com- 
pact presentation of these rates (i.e. for uo <C 2 ) for the three polarization modes 
of splitting permitted by CP (charge-parity) invariance in QED, namely _L— >-|||| , 
_L— >_L_L and II— >_L|| , is 



Jtf_„„(aO = B e sin% B Ml = ^5%, 
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using the superscript notation sp to signify photon splitting. Here the "scattering 
amplitudes," which have been computed numerically in [23l 124"] . are given by 

j> 1 f°° ds - S /B J ( 3 s\ coshs 3 + 2s 2 s cosh s\ 

1 ~ WJo T e \\ 4^ + 6 J ^ihhi + 12 sinh 2 s 2 sinh 3 s j ' 

(8) 

M 1 r°° ds £ - s /B f j coshs | 3 -As 2 3s 2 ] 

2 ~ B 4 Jo s I 4s sinh s 4 sinh 2 s 2 sinh 4 s \ 



Evaluations of these forms in terms of special functions were obtained in [31]: 




where T(x) is the Gamma function, ip(x) = dlog e T(x)/dx is the digamma function, 
L\ « 0.24875 is a constant resulting from the sum of the logarithmic series that 
expresses the remainder for Stirling's approximation for the T function, and 

log e ri(x) = jTdtlog e r(t) + ix(x-l)-|log e 27r . (10) 

At low fields, A4i ~ 26/315 and M. 2 ~ 48/315 are independent of B , but at high 
fields possess A4\ ~ l/(6£> 3 ) and M. 2 ~ l/(35 4 ) dependences. Hence, the rate for 
-L— HIH asymptotically approaches a constant as B — ► oo , a result that is exhibited 
in Figure [TJ These rates are of broad applicability to pulsar modeling. Deviations 
from this low energy limit near pair creation threshold are detailed by [231 H23 134"] : 
they are somewhat apparent in Figure [TJ for which the numerics were obtained 
from the full Schwinger-type results of [23] • Note that the B 6 dependence of the 
rates when B < 1 permits the use of so-called hexagon Feynman diagrams when 
treating the field as a small perturbation at three vertices, a technique employed in 
early computations of this process. This approximation cannot be used for higher 
B where the field's influence on the electron propagators must be incorporated 
self-consistently, and results in greater complexity as appears in Eq. (|8]). 

As mentioned in Section [2j the different pair creation attenuation rates for the 
two linear polarization states automatically generate polarization-dependent refrac- 
tive indices ( n± ^ nn ) and light speeds, i.e. birefringence. This birefringence of the 
magnetized vacuum also implies an alteration of the kinematics of strong field QED 
processes [23], admitting the possibility of non-collinear photon splitting. Hence, 
while the splitting modes _l_ — 1 1 , ||— >_L± and ||— are forbidden by CP invari- 
ance in the limit of zero dispersion, dispersive effects guarantee a small but non-zero 
probability for the _l_ — > _l_ 1 1 channel. Extensive discussions of linear dispersion in a 
magnetized vacuum [2U [35] demonstrate that in the limit of weak linear vacuum 
dispersion (roughly delineated by -Bsin^kB & B CT ), where the refractive indices for 
the polarization states are very close to unity, energy and momentum could be si- 
multaneously conserved only for the splitting mode -L— H| || (of the modes permitted 
by CP invariance) below pair production threshold. 



This result, known as Adler's [23] kinematic selection rules for photon splitting, is 
contingent upon the inequality R p f < R pp , assuming that pair creation 7 — > e + e~ 
provides the dominant contribution to the vacuum dispersion, generally true for 
B < 1 regimes of low plasma density]]] Therefore, it is probable that only the one 
mode (_L— >|||| ) of splitting operates in normal pulsars. However, this constraint 
may not hold in supercritical fields where stronger vacuum dispersion arises and 
other contributions may come into play. For example, the generalized vacuum po- 
larizability tensor may yield significant corrections from quadratic (i.e. those that 
connect to photon splitting) and higher order contributions. Recent analysis [37] of 
vacuum dispersion induced only by pair attenuation has indicated that while the 
selection rules are generally upheld as u increases towards pair threshold, they have 
to be modified in certain phase spaces above pair threshold. In such regimes, split- 
ting cannot be discounted relative to 7 — > e + e~ as a higher order process, since its 
rates sample pair resonances when the electron propagators go "on-shell," and the 
splitting rate may approach that of a first-order process. Such calculations of rates 
for 7 — > 77 at uj > 2/ sin^kB have not yet been performed at any length, though 
the formulations of [32j [33] can be readily applied to this problem. 

4 Anticipated Astrophysical Signatures 

Having encapsulated the essentials of the physics of these two exotic predictions of 
high-field quantum electrodynamics, the focus now turns to how signatures of these 
processes might realistically manifest themselves in pulsars. It is quite possible that 
a vindication in astrophysical contexts for the theoretical expectations may predate 
the eventual first observation of either process in terrestrial laboratories. 

4.1 Pair Creation 

Since the magnetic pair creation rate is generally a strongly increasing function of 
photon energy, its major signature is an attenuation turnover in pulsar radiation 
spectra. The resulting maximum energy of emission is controlled by attenuation 
during photon propagation through the pulsar magnetosphere. Such attenuation 
provides a characteristic super-exponential turnover [38] (see below) that contrasts 
that expected in outer gap models (e.g. see [SHI S3] for a comparison), scenarios 
where the gamma-rays are generated in a pulsar's outer magnetosphere and the 
quantum effects discussed here are negligible. Pair creation effectively occurs at 
the threshold cusin^kB = 2 for high fields, i.e. B ^ 1, and above threshold at 
u;sin#kB ~ 0.2/B for lower fields [15]. Hence, the mean free path for photon atten- 
uation in curved fields is A pp ~ p c /u> max{2, 0.2/B} , i.e. usually when usmO^B 
crosses above threshold during propagation. This assertion derives from the fact that 



1 Note that in the dense gases encountered deep in neutron star atmospheres, dispersion becomes 
dominated by plasma contributions [36J , for which the selection rules of Adler are not applicable. 



photons are generally emitted almost parallel to the local B due to the relativistic 
beaming associated with electrons/pairs streaming along the field lines. The radius 
of field curvature at altitude R > R^s is p c = [P Rq c/27r] 1//2 for a pulsar period 
P . This introduces a dependence of the energy when pair threshold is crossed of 
uj oc l/6*kB oc p c oc \/P~R~q . In a dipole field geometry, the local magnetic field scales 
as B oc Rq 3 , thereby generating the scaling uj oc 1 /B oc R$ for the energy of pho- 
tons that access pair threshold. Accordingly, the approximate dependence of pair 
creation cutoff energies e MAX on the surface polar field strength B , altitude Rq 
and pulsar period P (in seconds) can be summarized in the approximate empirical 
relation (e.g. jlO] and references therein) 

pf^V 72 max^ I. ( /,"" i ^> OV . (11) 
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Accurate numerical determinations of the turnover maximum energy from the codes 
developed in [UJ [T7j, are plotted in Figure [2] these include the effects of general 
relativity on spacetime curvature, field enhancement and photon energy in "slowly- 
rotating" systems with P ^> 2ttR ns /c. At fields B ^ 0.7 photon splitting acts to 
further reduce e MAX , as discussed in [T7j, though not shown in the Figure here. 
Figure [2] and Eq. (11) clearly indicate a strong anti-correlation between the max- 



imum energy and the surface magnetic field. Coupled with observed maximum en- 
ergies in pulsars detected by the Compton Gamma- Ray Observatory (CGRO), this 
behavior seems to be augmented by an apparent decline of emission altitude with 
B Q . Such a trend, which is not anticipated in outer gap models, is a distinctive char- 
acteristic that can be probed by NASA's upcoming Gamma-Ray Large Area Space 
Telescope (GLAST, to be launched in mid-2008). This core NASA astrophysics 
mission is expected to add at least 50-100 pulsars to the extant database, possibly 
many more, perhaps around half of which will offer clean e MAX determinations to 
refine Figure [2] further. The maximum energy is generally in the 1-10 GeV band for 
normal young pulsars, but can be much lower [1TI [T7| for highly magnetized ones, 
and also much higher for millisecond pulsars, so that signals in the 30-100 GeV 
band are possible [121 H3] for polar cap models via synchrotron/curvature cascades 
if the field is low enough. 

Such a statistical correlation for a population of pulsars is not the only signature 
of the operation of 7 — > e + e~ . Normally, maximum energy turnovers resulting from 
a moderately energy- dependent attenuation of photons in astronomical sources as- 
sume an exponential form <in 7 /cfe 7 oc exp{— e 7 /e cuto fr} • In this particular instance, 
the attenuation process is exponentially-sensitive [as evinced in Eq. Q] to the in- 
coming photon energy for subcritical field strengths; a similarly rapid onset should 
be experienced when B <; 1 . The upshot is that magnetic pair creation should 
impose a super- exponential cutoff in pulsar gamma-ray spectra of the approximate 
form dn^/dery oc exp{— a exp[— £ MAX /£ 7 ]} for some constant a . The severity of the 
cutoff will depend on the range of altitudes and magnetic latitudes sampled in the 



Pair Creation Cutoffs in Pulsars 
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Fig. 2: Maximum pulsar emission energies (adapted from [10]) imposed by pair 
creation attenuation at two different altitudes, Rq (dashed curves) and 3i?o (solid 
curves), described empirically via Eq. (If). For each altitude, a range of pulse 



periods (polar cap sizes) O.f — 5 sec is represented, as indicated. These energies are 
determined by the comprehensive photon propagation/attenuation code described in 
[17] . which includes curved spacetime effects. Observed/inferred cutoff energies (or 
ranges) for 8 gamma-ray pulsars of different Bq are indicated, from which a trend 
of declining altitude of emission with increasing surface polar field B is suggested. 



emission region, but the net result is typically [12l [3H] a sharp turnover that is dis- 
tinguishable from those generated by other attenuation processes. It is anticipated 
that GLAST will clearly be able to discern such super-exponential cutoffs, if they 
are present, in a number of pulsars, and accordingly offer discrimination between 
the viability of polar cap and outer gap models for high energy pulsar emission [44 J . 

The presence of strong fields virtually guarantees a strong polarization signal 
in polar cap models, and when these couple with spectral structure and temporal 
information, particularly powerful observational diagnostics are possible. Due to the 
polarization-dependence of the pair production rates, this may provide a potentially 
fruitful tool, particularly for highly-magnetized pulsars with B <^ f 13 Gauss, since 
the attenuation cutoffs fall in the 1-10 MeV band. Hard gamma-ray pair production 
tracking experiments like GLAST are generally not afforded the opportunity to act 
as polarimeters, being limited by multiple scattering in tracking chambers above 



300 MeV. Medium energy gamma-ray experiments, on the other hand, are ideally 
suited to polarization studies, essentially via their sampling of Compton scattering 
kinematics. Gamma-ray polarimetry is no longer a distant dream [33]. Europe's 
current INTEGRAL mission, and NASA's solar RHESSI telescope have some ability 
to detect polarization in bright sources such as gamma-ray bursts (IHlHZj and X-class 
solar flares. Moreover, polarimetric capability in the hard X-ray and soft gamma-ray 
bands is a high priority for next-generation tracking Compton detectors such as the 
Advanced Compton Telescope experimental initiative [IS]. 

As an interesting aside pertaining to a completely different astrophysical context, 
it is worth noting that opacities for pair creation are currently included in simulation 
codes for the propagation of ultra-high energy cosmic rays and air shower initiation 
in the Earth's upper atmosphere and magnetosphere. This presumes a possibility 
that these energetic particles could be photons, something that is currently not a 
popular perception, but has not yet conclusively been excluded. For completeness, 
the shower simulations for ground-based detector arrays (e.g. [IHJEO]) such as Fly's 
Eye/HIRES, AGASA and Auger include 7 — > e + e~ rates when treating photon 
primaries, which for terrestrial fields of B ~ 1 — 10 Gauss sample the asymptotic 
regime in Eq. ([6|. Any photons in the e 7 ~ 10 20 eV domain can effectively convert 
into pairs off the geomagnetic field above the ionosphere and thereby instigate a 
so-called "pre-shower" that is subsequently reprocessed in a complicated hadronic 
cascade in the atmosphere below. Such a convolution of interactions precludes defini- 
tive determination of the action of 7 — > e + e~ even if a photon primary origin were 
eventually concluded. Yet a remarkable aspect of this field is that the participating 
physicists attribute much greater certainty to the strong-field QED physics than to 
extrapolations of B = photonuclear cross sections to such extreme energies. 

4.2 Photon Splitting 

It is clear from the rates displayed in Figure [l] that the necessary requirement for 
photon splitting to operate in pulsar magnetospheres is that B ^ \ . This turns the 
focus to a small subset of known neutron stars/pulsars that possess ultrastrong fields, 
termed magnetars. Observations at X-ray energies have yielded detections of both 
long periods and high period derivatives in two types of such sources, anomalous 
X-ray pulsars (AXPs) and soft 7-ray repeaters (SGRs), which suggest dipole spin- 
down fields in the range 10 14 — 10 15 Gauss. The AXPs are a group of seven or eight 
pulsating X-ray sources with periods in the range 6-12 seconds, and are continuously 
spinning down [51J. The SGRs are a type of 7-ray transient source that undergoes 
repeated bursts, some of them enormous (e.g. [52]); four are definitively known to 
exist. With the possible exception of the unconfirmed reports [53] of a low frequency 
pulsed detection of a counterpart (PSR J1907+0919 at 111 MHz) to the soft gamma 
repeater SGR 1900+14, none of these sources has detectable pulsed radio emission. 
It is notable, however, that an unusual related source, the transient AXP XTE 
J1810-197, has exhibited transient and peculiar pulsed radio emission [54] . 



A central principal of the radio pulsar paradigm is that the emission is coupled 
to a prolific presence of pairs [HE]. As pulsars age, their periods lengthen, due 
to the action of magnetic dipole radiation torques. Concomitantly, the geometrical 
confines for their acceleration and emission locales on open field lines become more 
constricted, roughly as the polar cap opening angle 9 cap fa [27tPns/(-Pc)] 1//2 shrinks. 
This establishes a beaming of radiation generated by accelerated pairs, and the 
collimation angles #kB for potentially pair producing photons drop. Eventually, the 
traversal distance in the magnetosphere required to cross pair threshold exceeds the 
scale for magnetic field decline (i.e. ~ Pns)- P arr creation then shuts off, and 
presumably so does the radio emission, i.e. radio pulsars "die." This has led to 
the concept of a pulsar death line at long periods, something that is illustrated in 
Figure^ which is the P-P diagram, the conventional depiction of pulsar phase 
space. The location of this boundary depends on observables P and the period 
time derivative P , through the relation for the polar field B Q that is inferred for 
magnetic dipole torques seeding rotational deceleration (e.g. [T7]): 

B Q = 6.4 x 10 19 [PP] 1/2 Gauss . (12) 

Dotted diagonal lines in Figure [3] denote contours of constant field strength. The 
locations of over 1500 observed radio pulsars are marked in this Figure. Three of 
the 8 known gamma-ray pulsars, namely the Crab, Vela and PSR 1509-58 are also 
highlighted (see lower inset), together with the positions of five radio-quiet anoma- 
lous X-ray pulsars (filled green triangles) and SGRs 0526-66, 1806-20 and 1900+14 
(filled red squares) in the upper right, constituting the sub-class of magnetars. 

An alternative mechanism for suppression of radio emission, specifically in high 
field pulsars and magnetars, was proposed by [551 Ej- F° r local fields B <; 1 , pho- 
ton splitting can be an effective competitor to pair creation. This arises principally 
because the pair threshold is crossed only after emitted photons are transported 
through the magnetosphere to acquire significant angles with respect to field lines. 
During this time, when B ^> 1 , photon splitting can have time to act, and its prob- 
ability of attenuating the photons and pushing them further below pair threshold 
by energy degradation is sensitive to the local field strength B and the polar cap 
geometry, i.e 6 cap . Consequently, the competition between splittings _L— >|||| and 
pair conversions e + e~ is defined by the principal pulsar observables P and P . 
Detailed discussions and computations of this competition are presented in [35] [TT]. 
They include all the features of QED discussed above, and the general relativistic 
distortions of the magnetosphere and the photon trajectories in a Schwarzschild ge- 
ometry. The principal result is an identification of the approximate P and P pulsar 
phase space when it is anticipated that photon splitting will dominate pair creation 
for _L photons. This is depicted in Figure |3j Since pair suppression is putatively 
identified with a lack of radio emission in pulsars, this defines regimes where photon 
splitting might seed radio quiescence in SGRs and AXPs. Clearly there are some 
recently-discovered pulsars above the boundary, a situation contrasting the complete 
demarcation highlighted in [55] that pertained to the Princeton Catalogue era. 




Log 10 [Period, P (sec)] 



Fig. 3: The upper portion of the P-P diagram, with filled circles denoting the 
locations of 541 members of the Princeton Pulsar Catalogue [56], and open circles 
marking 1599 radio pulsars in the more recent ATNF Pulsar Catalogue [57] (only 
those with P > 10 -15 are visible). The approximate boundary (diagonal red line) 
demarcating when splitting dominates (uppermost P ) and when pair creation is 
prolific (lower P for conventional pulsars) is indicated, being computed in |17j . 
Above this, the heavy-weight curves labelled e esc = 10 MeV and e esc = 10 2 MeV are 

contours for the escape energy of photon splitting, I — ^ ] 1 1| ; to the right of these, the 

magnetosphere is transparent to 10 MeV and 100 MeV photons, respectively. 

The placement of the quiescence boundary in the Figure presumes that the emis- 
sion region is very near the stellar surface, a criterion that was specifically adopted 
when identifying the boundaries in Figure [3} If the altitude Ro of the photon atten- 
uation locale moves above the surface, then the putative radio quiescence boundary 
correspondingly moves to higher P , i.e. higher surface polar fields B . Since the 
local dipole field scales as B ~ B (R /Rns)~ 3 , which is pinned by the balanced 
competition between pair creation and splitting, and the inferred value of B scales 
with observables as Bo oc (PP) 1 ' 2 , the location of the quiescence boundary in the 
P-P diagram approximately obeys P oc R$ 6 for fixed pulsar periods P (and there- 
fore also fixed polar cap angular sizes). Hence, small increases in R above i?NS can 
easily [UJ move the boundary above the radio pulsar population displayed in the 



Figure. It should be emphasized that if Adler's splitting selection rules [23] apply 
here, then photons of || polarization are only subject to pair creation. This implies 
only a partial suppression of 7 — > e + e~ by splitting, which may lead to a post- 
ponement of pair creation eventually to higher altitudes [T7j. However, photons of 
_L polarization dominate the photon population generated by the principal primary 
and secondary emission processes of polar cap models for 7-ray pulsars [12, 38J, 
namely cyclotron/synchrotron radiation, curvature emission and resonant Comp- 
ton scattering. Hence, the partial suppression of pair creation by photon splittings 
_L— is in fact likely to be quite effective at typical magnetar field strengths. 

While this identifies a signature and observable action of photon splitting in 
an astrophysical setting, the connection to radio emission is the weak link in this 
logical construct. The commonly- held belief that the presence of a profusion of 
pairs is requisite for strong radio emission is not yet proven. The actual nature 
of the coherent radio emission process(es) is still unknown, a major problem in 
the study of pulsars. If the density of pairs is not relevant for the production of 
bright radio emission, then the action of splitting in magnetar field regimes has no 
formal connection to the observed radio quiescence or faintness of SGRs and AXPs. 
Hence this property cannot, at present, be regarded as a definitive measure of the 
action of splitting in pulsars. For potentially more direct evidence for the operation 
of 7 — > 77 in pulsars, one turns again to spectroscopic evidence similar to that 
discussed in Section |4.1| for pair creation. Noting the maximum energy contours in 
Figure[3j for magnetars and pulsars with fields above 4x 10 13 Gauss, photon splitting 
should prohibit any emission above ~ 100 MeV, though prominent signals below 
100 MeV are possible [5B] due to the efficiency of resonant Compton upscattering of 
surface thermal X-ray photons. Accordingly, photon splitting should generate the 
spectral turnovers in supercritical fields, not solely pair creation, and the physics 
of the process guarantees that the cutoff energies are polarization-dependent. This 
effect was discussed in |JT], where a focal study of the high-field pulsar PSR 1509-58 
indicated that the observed spectral data in the sub- 100 MeV range implied a cutoff 
that could only be modelled in a polar cap scenario if photon splitting was acting in 
addition to general relativistic effects. The data for this pulsar do not permit much 
play in the parameter space, so the inference that splitting is active in PSR 1509-58 
is quite strong. Yet it is crucial to find one or more other high- B pulsars where such 
constrained inferences can be made. Unfortunately such a prospect is not yet on 
the horizon. While NASA's imminent GLAST mission will discover an array of new 
gamma-ray pulsars, this issue can only be addressed by a sensitive, next-generation 
telescope in the 1-30 MeV range such as the Advanced Compton Telescope. Yet 
this specific science question, acquiring strong evidence for the action of an untested 
prediction of strong-field QED, forms part of the many motivations for developing 
a new medium-energy gamma-ray mission that has polarimetric capability. 



5 Conclusion 



Two of the exotic and fascinating predictions of quantum electrodynamics in strong 
magnetic fields have been highlighted in this paper, one-photon pair creation and 
photon splitting. It is an impressive indicator of the confidence placed in the theory 
of QED by physicists that, in spite of a lack of experimental vindication in terrestrial 
laboratories, they are both calculated theoretically in spite of significant degrees of 
mathematical complexity (particularly for photon splitting), and anticipated almost 
routinely in the astrophysical context of pulsars and elsewhere. There have been 
various limited attempts to explore these processes in laboratories, sometimes indi- 
rectly by searching for birefringence signatures using lasers in optical cavities [59, 60J, 
but so far without definitive results. Largely this is because the required fields are 
beyond present experimental capabilities: for B ~ 1 MGauss, super- TeV energy 
photons are required to start to approach the criterion ujB ^0.1 (in dimensionless 
units) for 7 — > e + e~ to become significant. It is not beyond reason that building 
from the status quo of current powerful terrestrial magnets and particle accelerator 
technology such as the Large Hadron Collider (LHC), it may prove possible to access 
the appropriate u - B phase space in a laboratory setting in the coming decades. In 
the meantime, it is of great interest to explore their application to astrophysics, not 
just to anticipate an alternative proof of the action of these QED processes, but also 
to enhance our understanding of the exotic cosmic sources, specifically pulsars, that 
can provide a conducive environment for the activity of 7 — > e + e~ and 7 — > 77 . 
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